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We study a double quantum dot in the regime where each dot carries a spin-1/2. This system 
is described by the 2-impurity Kondo model, having a non-Fermi liquid fixed point for a critical 
value of the inter-impurity coupling. The Hamiltonian describing the vicinity of the critical point, 
including the relevant potential scattering perturbations, can be cast in quadratic form. This allows 
us to predict a universal scaling function for the finite temperature nonlinear conductance along the 
crossover from the critical point to the surrounding stable fixed points. 
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Introduction. It is now well established that quantum 
dots (QDs) behave as Kondo impurities at low temper- 
atures [l|, [2|. So far, most experiments focused on the 
study of Fermi-liquid states, with regular thermodynamic 
as well as nonequilibrium transport properties 3 at low 
temperatures. However, recently a controlled QD setup 
was constructed experimentally to access the two-channel 
Kondo fixed point 4]. Remarkably, the conductance of 
this fixed point matches a non-Fermi-liquid (NFL) scal- 
ing function over a broad range of energy. 

Another simple impurity model showing NFL behav- 
ior, similar to the two-channel Kondo model, is the two- 
impurity Kondo model (2IKM), consisting of two im- 
purity spins that are coupled to conduction electrons 
and, at the same time, interact with each other through 
an exchange interaction K. Jones et. al. [f| observed 
that in the 2IKM a NFL quantum critical point (QCP) 
at K = K c separates a local singlet from a Kondo- 
screened phase. The exact critical behavior was found 
using conformal field theory methods in @, 0| • Georges 
and Meir [§] and Zarand et. al. [§] proposed different QD 
realizations of this critical point. 

In this paper we study a series double quantum dot re- 
alizing the 2IKM when each dot behaves as an effective 
spin [see Fig. ©]. Ref. 1 studied this system in a wide 
range of parameters at zero temperature using a slave bo- 
son mean field theory. Here we concentrate on the vicin- 
ity of the critical point and predict exact full crossover 
functions for the nonlinear conductance from the QCP 
to the surrounding manifold of fixed points at finite tem- 
perature and source drain voltage. This crossover is cal- 
culated taking into account all three relevant perturba- 
tions at the QCP. The first perturbation is associated 
with K — K c . Whereas typically potential scattering is 
a marginal perturbation in quantum impurity models, in 
the 2IKM it leads to two additional relevant operators at 
the critical point @, 0, H| ■ 

Our results are obtained using a convenient basis in 
which the QCP is described by a quadratic Hamilto- 
nian. This Hamiltonian follows from Gan's theory for 
the 2IKM [1(3]. While this theory was derived for a spin 



anisotropic version of the 2IKM, it was argued that it 
describes correctly also the spin SU(2) symmetric criti- 
cal point [l(| ■ We substantiate this statement further by 
showing explicitly that the operators at the critical point 
have the same form as in the conformal field theory . 

In our crossover calculation of the conductance we in- 
clude not only finite temperature, but also address the 
non-equilibrium problem at finite source drain voltage. 
Originally this problem was addressed using perturba- 
tion theory [111 ], valid at energy scales large compared 
to the Kondo temperature Tk- To address the crossover 
to low energies different nonperturbative techniques were 
adopted. Ref. [L2j studied transport through a 1-channel 



Kondo impurity using abelian bosonization; however ex- 
act results were obtained only for a specific point in 
the parameters space (Toulouse limit). Another impor- 
tant development in this direction was the application 
of the Bethe-ansatz and finding of many-body scattering 
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FIG. 1: Double quantum dot system. 



Model. Double QDs have already been proven to dis- 
play rich Kondo physics in competition with inter-dot ex- 
change interaction [1, IE EE EE E3, El ■ A multilevel sin- 
gle QD near a singlet triplet transition is also described 
by the 2IKM 

Our model consists of two leads, denoted L and R (or 
i = 1 or 2 respectively), attached to two QDs which 
behave as effective spins Sr,, Sr due to large Coulomb 
energy U. After the standard "unfolding transforma- 
tion" [2(| , one obtains a description of the model in terms 
of 4 left moving Dirac fcrmions ipi a (x), i — L, R, a =t, |, 
x 6 {— co, oo}. We use a convention in which H = Vp = 1, 
where vf is the Fermi velocity. At the moment let us 
assume particle-hole symmetry. The interaction Hamil- 
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tonian generated to leading order in the tunneling am- 
plitudes ii, tn, takes the form 



Hr — Jlsll ■ Sl + Jrsrr ■ Sr + KSl • Si- 



(1) 



with Sij — ip^i^f-ipjp, Ji ~ tf/U, and K ~ t 2 LR /U. 

At tifl = the Kondo effect takes place separately 
between the left lead and left dot, and between the 
right lead and right dot. The system is insulating in 
this case. This Kondo-screened phase persists at pos- 
itive nonzero K up to K = K c , above which a local 
singlet forms between the two impurities [H, 0|. The 
critical value of K occurs at K c ~ Tk- Here we as- 
sume roughly equal Kondo temperatures for left and right 
sides, Tkl ~ Tkr = Tk, but most of our results do not 
depend on that assumption, as discussed below. The 
fixed point at K = K c has NFL properties and is unsta- 
ble with respect to an operator with anomalous dimen- 
sion A = 1/2 associated with the perturbation K — K c . 
However for the general case when particle-hole and par- 
ity (L <-> R) symmetries are broken the QCP is unstable 
against two additional relevant operators identified with 
potential scattering @, 0] • It was found that while intra 
lead potential scattering terms Vll^liPl + VrriPr^r are 
marginal at the QCP, inter-lead direct tunneling terms 
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V L Rlp L 1pR + h.c. 



(2) 



correspond to two relevant A — 1/2 operators. In our 
model \Vlr\ ~ titLR.tR/U 2 . When parity is broken Vlr 
can be complex. 

Fixed point Hamiltonian. We bosonize the original 
theory and introduce four left moving bosonic fields <pj a 



satisfying ipj a 



oi<t>i 



Subsequently 4 linear bosonic 



combinations are defined, corresponding to charge, spin, 
flavor, and difference of spin between the flavors: 4> c — 

lEja&ai <f>s = |Eja(°" Z )S0jc«. <t>f = \ Y.jJf)j<l>i>*> 

and (f) X = |Lj Q ( r2 )j( cr2 )a^a- Here & and ? are 
Pauli matrices acting on the spin and L/R space re- 
spectively. Since the exponentials of the new bosons 
have dimension A = 1/2, we can define new fermions 



— ^9 A A — 



, A = c,s, /, X . Taking the real and imaginary 



X2 



o nly 

X^(0+) [I3,l2l|. This follows by identifying 



it is found that the o nly nontrivial part of the BC reads 
x with the Majorana fermion of the Ising model in 



It is our main point in this section to argue that for 
energy scales well below Tk the corrections to the free 
Hamiltonian of the QCP take a simple quadratic form in 
the 5*0(8) representation: 

3 

SH = iJ2*jXj(0)a, (3) 

K — K c 

Ai =ci — 7=A (A 2 ,A 3 ) = c 2 \/TKv(ReV L R,lmV L R), 



where Xi( x ) = X? (^)sgn(x), {x2,Xs) = [xLx{)> c i and 
C2 are constant factors of order 1, and a is a local Majo- 
rana fermion a 2 = 1/2. 

At A2 = A3 = 0, Eq. corresponds to the correction 
to the fixed point hamiltonian of the particle-hole and 
parity symmetric model, Eq. ([1]). Using conformal field 
theory methods it was shown that the relevant operator 
in this case is equivalent to a magnetic field acting on the 
boundary spin of a quantum Ising chain. This is consis- 
tent with the j = 1 term of E q. (|3l ) , since one can identify 
Xi(x) with the Ising fermion [21] and write the boundary 
spin of the Ising model as ob — iXi( x = 0) a "SI]- I n this 
analogy Ai corresponds to the magnetic field acting on 
the boundary spin. 

In addition we will clarify that the j = 1 term of 
Eq. ([3]) is consistent with Gan's theory for the anisotropic 
Toulouse limit of the 2IKM [nj]. This theory uses the 
50(8) Majorana representation, and the two impurity 
spins turn into a local fermion d, where {d, d^} = 1. 
Defining two Majorana fermions a = and b = , 
the boundary part of Gan's Hamiltonian reads 



SHr 



:X2 (0)6 - i{K - K c )ab. 



(4) 



In general SHq 7^ SH\\ 2= \ 3= q, however we shall show 
that the two coincide for energy scales <§; Tk- To see 
this suppose K — K c and consider a mode expansion 



parts of those fermions we obtain 8 Majorana fermions X2 ( x ) = / [Vki&tyk + h.c ), b = / (ukipk + h.c), 



A _ ^ +i>A A 
Xl ~ 7f , X2 



l/l T — 1pA 

V2i 



This is denoted as the 



50(8) Majorana representation of the free fermions, due 
to the 50(8) symmetry group at the trivial fixed point 
J = 0. 

At the trivial fixed point J = the free Majorana 
fermions have trivial boundary condition (BC) x^(0~) = 
xt(® + ), i = 1, 2, A = c, s, /, X, implying the continuity 
of those fields at the boundary x = 0. The fermions are 
still free at the critical point and obey the BC estab- 
lished in Ref. [6(, using a Bose-Ising representation for 
the free fermions. It was found that the nontrivial BC 
occurs in the Ising sector of the theory. By relating this 
representation with the 50(8) Majorana representation, 



where {ipk, "0fc<} — S(k—k'), {tpk, ^fe'} = 0, and where ini- 
tially we choose A 3> Tk as an ultraviolet cutoff. Solving 
Schrodinger's equation for the wave functions <fk(x) and 



Uk one finds tpk {x) 



<^fc(0) 



„2iS 



tan S = One should normalize = where £ 

is the size of the system. While at Tk = we have the 
BC X2"(0 + ) = Xa we see from the wave function 

that the effect of the first term in SHq is to modify this 
BC to X2 (0 + ) = -X2 (0~) for energies < T K . The key 
observation is that the following operator identity holds 



3 



if one restricts the mode expansion of its LHS and RHS 
to energies below a cutoff A <C Tk, 

1 



b = 



(5) 



Physically this means that at energy scales below Tk the 
local operator b is absorbed into the field X2 an d changes 
its BC. Using the operator identity Eq. (JSJ), we see that 
the term oc K — K c in 5Hq is equivalent to the j = 1 
term in 5H. This establishes the connection between 
Gan's theory and the boundary Ising model arising from 
the conformal field theory solution, showing that Gan's 
anisotropic theory describes correctly also the vicinity of 
the isotropic fixed point. 

When particle-hole and parity symmetries are broken 
Gan's theory predicts two additional relevant A = 1/2 
operators which take the form of the j = 2,3 terms 
of Eq. ©. The presence of these two terms is con- 
sistent with the conformal field theory: Consider the 
relevant potential scattering terms ip L ipR + ipRipL an d 
i{^\^R ~ tpRipL) in Eq. which have the 50(8) rep- 
resentation 1X2 Xi an d 1X2 X2 respectively. We are in- 
terested in the form these operators take near the NFL 
critical point. This can be obtained, using conformal field 
theory methods, by applying double fusion in the Ising 
sector. This operation acts on \2 ov ^li an( A transforms 
it into the identity operator. To maintain the bosonic 
nature of the Hamiltonian a local dimension Majorana 
fermion must replace \2 ■ To account for the correct 
ground state degeneracy of the QCP one requires this lo- 
cal Majorana fermion to coincide with the same operator 
a which couples to X2 K ^ K c . Therefore the two 
relevant potential scattering terms must have the form of 
the j = 2, 3 terms of Hqcp- 

Conductance. With the fixed point Hamiltonian 
Eq. ((3]) we proceed to calculate the conductance of the 
double dot along the crossover from the critical point to 
the surrounding fixed points, as function of temperature 
T, source drain voltage V, and the three coupling con- 
stants Ai, A2, and A3. As we shall see, the crossover is 
controlled by an energy scale 



T* = A 2 + A 2 + A 2 



(6) 



Our result, Eq. ©, is valid for any T/T* and eV/T*, as 
long as eV,T, T* < min{7W,, T K r}. 

Schiller and Hershfield [L2J studied the related prob- 
lem of transport through a single QD with effective spin 
5. They considered the anisotropic Kondo interaction 
of the form Hjf = T,i,j=L,R 4* a l SV with 

J l J — Jy ] — Jj/. Using the 50(8) representation, they 
obtained a quadratic theory [Eq. (3.15) at zero magnetic 
field B = 0] for the case J 2 ii? = 0, J z Li = jf R = 2tt, 
with a boundary term 

:[J + X 2 xb + J^ R X 1 f a + J-j&a], (7) 



where = (Jf L ± J RR )/2 and a -1 is an ultraviolet 
momentum cutoff. Comparing 5 Hid to 5Hq in Eq. 
the coupling J + Xx^ * s t Qe only term involving b. For 
either the double or single QD the current operator is 
given by I = -i[H, Y] where Y = \ di(f[fe - 

V ' RtpRa) — f^odxi/jfipf. We see that the current op- 
erator involves the flavor fermions, and hence the first 
term oc J + in Eq. does not contribute to the con- 
ductance, since the flavor fermions are decoupled from 
b. The operator a in Eq. couples linearly to bulk 
Majorana fields at x = 0. Comparing to Eq. ([3]), it is 
easy to see that under the replacements £ 7ro Ai and 

-^== <-> y \\ + A3 , the two systems have the same cur- 
rent as function of temperature and voltage. We shall 
not repeat the calculation of [12| and present the result 
for the nonlinear conductance, G = dl/dV, in terms of 
the parameters of our model 



G — GqF 



T 



eV 
1 



Go 



2e 2 T 



LR 



h T* 



F[t,v] = -Re^' 



1 1 iv \ 

2 + 2t + 2^t) 



(8) 



where T£ H = A| + A3 and ip{ z ) is the digamma function. 
We see that the energy scale T* defined in Eq. ([6]) deter- 
mines the crossover scale. At K — K c , T* — > T£ fi . The 
fact that T* is quadratic in K — K c and in Vlr is a signa- 
ture of a NFL; these coupling constants have fractional 
renormalization group scaling dimensions 1/2. In Fig. @ 
we plot the conductance as function of eV/T* for various 
values of T/T* . At T = 0, G/G = (l + (eU/TrT*) 2 ) -1 
and for V = 0,T< T* , G/G -> 1 -T 2 /3T* 2 . 




eV/T* 

FIG. 2: Scaling function for the nonlinear conductance. 

Note that the current vanishes when Vlr = 0. This 
occurs since the corresponding operators in Eq. @ are 
the only terms in our model which do not conserve Y. At 
the fixed point we expect additional terms not conserving 
Y, such as the marginal operators X12X2 ■ However the 
amplitude for such terms transferring charge between the 
leads is small, and their effect can be neglected. 

Experimentally it is convenient to tune K through 
K c by varying tLR, and observe the evolution of a con- 
ductance peak as function of energy scales T and eV. 



4 



At K = K c there is a residual energy scale T£ R due 
to potential scattering. Using \Vlr\ ~ tLt [}i tR and 
t 2 

if ~ ~ Tk we see that for a typical QD the en- 
ergy scale T^r is in the range of validity of our the- 
ory Tl R ~ (vJ) 2 T 2 /U « T K @. At low energies T, 
< T*, we have F — > 1 and G — > G . Thus G exhibits 
a peak in the low T linear conductance with height 2e 2 jh 
at K = K c , and width in K — K c which scales with Vlr- 
Such a peak structure was discussed in [8[; here we give 
its precise form. Intra-lead potential scattering Vll an d 
Vrr are marginal perturbations; they reduce the con- 
ductance at the peak to G pea k = — 0(1/ J) 2 ]. Now 
consider the regime T, eV ^> T*. In this case as t^R is 
varied the system passes through the critical region at 
K ~ if c where the system does not flow away from the 
NFL fixed point due to any of the relevant perturbations. 
The conductance now becomes: 



G 



e 2 T* 



Re V' 



ieV\ 
2tvTJ 



h T 2cosh 2 §^' 



(9) 



The width of the peak in G as a function of source-drain 
voltage, V, now scales with T. 

Although we assumed Tkl ~ Tkr above, in fact, 
most of our results don't depend on that assumption 
since we have already taken into account all relevant and 
marginal operators allowed in the effective Hamiltonian 
when parity is broken. Only the precise values of T*, 
T£ H and K c change in the asymmetric case, and the nec- 
essary conditions to approach the critical region becomes 
T* , T, eV -c min{TKR, Tkl}- In particular, note that 
the T = linear conductance has the value 2e 2 /h, for 
K = K c , independent of asymmetry in that limit. This 
is in striking contrast to the T = linear conductance 
through a single Kondo impurity which is suppressed by 
a factor of 2\tLtR\/ {\tL\ 2 + |tij| 2 ). In the strongly asym- 
metric case with Tkr <C Tkl, we find, up to logarithmic 
corrections, that 



T. 



LR 



T* « [K - K c ) 2 T K r/T kl + \vVlr\ 2 T K r, 
« \vV L r\ 2 T K r, and K c « T KL @]. 



(10) 



As pointed out by Schiller and Hershfield [12j . their 
single quantum dot model has the peculiarity of be- 
ing related to the two-channel Kondo NFL fixed point. 
This may be largely a consequence of fine tuning to the 
Toulouse limit J 2 Li? = 0, L = J* R = 2vr. The two- 
channel and the two-impurity NFL critical points are 
closely related. The double dot system discussed here 
provides a method of getting this NFL behavior in an 
experimentally feasible way rather than by the unrealis- 
tic fine tuning of parameters. 

Our system differs from the proposal of Zarand et. 
al. @, in the way in which the leads couple to the inter- 
acting system. In only one linear combination of the 
source and drain leads acts as a Kondo-screening chan- 
nel, whereas the second combination acts like a scan- 
ning tunneling microscope tip. A third lead is needed as 



a second Kondo screening channel. On the other hand 
the standard double dot system considered here has only 
two leads and can be therefore easier to construct ex- 
perimentally. In 0] exotic temperature dependence for 
the conductance \[T arises from the irrelevant operator. 
However we find that the (same) irrelevant operator gives 
regular quadratic corrections in eV and in T for our sys- 
tem [24j. On the other hand the NFL signatures are ap- 
parent from the relevant operators with dimension 1/2 
which can be addressed analytically in our configuration. 

Conclusion. In this paper we found new exact results 
for the 2IKM in a double QD. We derived full crossover 
formulas for the conductance from the NFL critical point 
to the stable fixed points in a 3- dimensional parame- 
ter space including K — K c and the two additional rel- 
evant potential scattering operators. We clarified that 
this crossover is described by a quadratic Hamiltonian. 
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